Programming Molecular Devices using Nucleic Acid Hairpins

Nucleic Acid hairpins have been a subject of study for the last four decades. They are composed of single strand that is

hybridized to itself, and the central section forming an unhybridized loop. In nature, they stabilize single stranded RNA, serve as nucleation

sites for RNA folding, protein recognition signals, mRNA localization and regulation of mRNA degradation. On the other hand,

DNA hairpins in biological contexts have been studied with respect to forming cruciform structures that can regulate gene expression.

The use of DNA hairpins as fuel for synthetic molecular devices, including locomotion, was proposed and experimental demonstrated in 2003. They

were interesting because they bring to the table an on-demand energy/information supply mechanism.

The energy/information is hidden (from hybridization) in the hairpin’s loop, until required.

The energy/information is harnessed by opening the stem region, and exposing the single stranded loop section.

The loop region is now free for possible hybridization and help move the system into a thermodynamically favourable state.

The hidden energy and information coupled with

programmability provides another functionality, of selectively choosing what reactions to hide and

what reactions to allow to proceed, that helps develop a topological sequence of events.

Hairpins have been utilized as a source of fuel for many different DNA devices. In this thesis, we program four different

molecular devices using DNA hairpins, and experimentally validate them in the

laboratory. 1) The first device: A

novel enzyme-free autocatalytic self-replicating system composed entirely of DNA that operates isothermally. 2) The second

device: Time-Responsive Circuits using DNA have two properties: a) asynchronous: the final output is always correct

regardless of differences in the arrival time of different inputs.

b) renewable circuits which can be used multiple times without major degradation of the gate motifs

(so if the inputs change over time, the DNA-based circuit can re-compute the output correctly based on the new inputs).

3) The third device: Activatable tiles are a theoretical extension to the Tile assembly model that enhances

its robustness by protecting the sticky sides of tiles until a tile is partially incorporated into a growing assembly.

4) The fourth device: Controlled Amplification of DNA catalytic system: a device such that the amplification

of the system does not run uncontrollably until the system runs out of fuel, but instead achieves a finite

amount of gain.

Nucleic acid circuits with the ability

to perform complex logic operations have many potential practical applications, for example the ability to achieve point of care diagnostics.

We discuss the designs of our DNA Hairpin molecular devices, the results we have obtained, and the challenges we have overcome

to make these truly functional.

Gene expression disruptions of organism versus organ in Drosophila species hybrids.

Hybrid dysfunctions, such as sterility, may result in part from disruptions in the regulation of gene expression. Studies of hybrids within the Drosophila simulans clade have reported genes expressed above or below the expression observed in their parent species, and such misexpression is associated with male sterility in multigenerational backcross hybrids. However, these studies often examined whole bodies rather than testes or had limited replication using less-sensitive but global techniques. Here, we use a new RNA isolation technique to re-examine hybrid gene expression disruptions in both testes and whole bodies from single Drosophila males by real-time quantitative RT-PCR. We find two early-spermatogenesis transcripts are underexpressed in hybrid whole-bodies but not in assays of testes alone, while two late-spermatogenesis transcripts seem to be underexpressed in both whole-bodies and testes alone. Although the number of transcripts surveyed is limited, these results provide some support for a previous hypothesis that the spermatogenesis pathway in these sterile hybrids may be disrupted sometime after the expression of the early meiotic arrest genes.

Morphological and genomic characterization of Filobasidiella depauperata: a homothallic sibling species of the pathogenic cryptococcus species complex.

The fungal species Cryptococcus neoformans and Cryptococcus gattii cause respiratory and neurological disease in animals and humans following inhalation of basidiospores or desiccated yeast cells from the environment. Sexual reproduction in C. neoformans and C. gattii is controlled by a bipolar system in which a single mating type locus (MAT) specifies compatibility. These two species are dimorphic, growing as yeast in the asexual stage, and producing hyphae, basidia, and basidiospores during the sexual stage. In contrast, Filobasidiella depauperata, one of the closest related species, grows exclusively as hyphae and it is found in association with decaying insects. Examination of two available strains of F. depauperata showed that the life cycle of this fungal species shares features associated with the unisexual or same-sex mating cycle in C. neoformans. Therefore, F. depauperata may represent a homothallic and possibly an obligately sexual fungal species. RAPD genotyping of 39 randomly isolated progeny from isolate CBS7855 revealed a new genotype pattern in one of the isolated basidiospores progeny, therefore suggesting that the homothallic cycle in F. depauperata could lead to the emergence of new genotypes. Phylogenetic analyses of genes linked to MAT in C. neoformans indicated that two of these genes in F. depauperata, MYO2 and STE20, appear to form a monophyletic clade with the MATa alleles of C. neoformans and C. gattii, and thus these genes may have been recruited to the MAT locus before F. depauperata diverged. Furthermore, the ancestral MATa locus may have undergone accelerated evolution prior to the divergence of the pathogenic Cryptococcus species since several of the genes linked to the MATa locus appear to have a higher number of changes and substitutions than their MATalpha counterparts. Synteny analyses between C. neoformans and F. depauperata showed that genomic regions on other chromosomes displayed conserved gene order. In contrast, the genes linked to the MAT locus of C. neoformans showed a higher number of chromosomal translocations in the genome of F. depauperata. We therefore propose that chromosomal rearrangements appear to be a major force driving speciation and sexual divergence in these closely related pathogenic and saprobic species.

Chromosomal organization of adrenergic receptor genes.

The adrenergic receptors (ARs) (subtypes alpha 1, alpha 2, beta 1, and beta 2) are a prototypic family of guanine nucleotide binding regulatory protein-coupled receptors that mediate the physiological effects of the hormone epinephrine and the neurotransmitter norepinephrine. We have previously assigned the genes for beta 2- and alpha 2-AR to human chromosomes 5 and 10, respectively. By Southern analysis of somatic cell hybrids and in situ chromosomal hybridization, we have now mapped the alpha 1-AR gene to chromosome 5q32----q34, the same position as beta 2-AR, and the beta 1-AR gene to chromosome 10q24----q26, the region where alpha 2-AR is located. In mouse, both alpha 2- and beta 1-AR genes were assigned to chromosome 19, and the alpha 1-AR locus was localized to chromosome 11. Pulsed field gel electrophoresis has shown that the alpha 1- and beta 2-AR genes in humans are within 300 kilobases (kb) and the distance between the alpha 2- and beta 1-AR genes is less than 225 kb. The proximity of these two pairs of AR genes and the sequence similarity that exists among all the ARs strongly suggest that they are evolutionarily related. Moreover, they likely arose from a common ancestral receptor gene and subsequently diverged through gene duplication and chromosomal duplication to perform their distinctive roles in mediating the physiological effects of catecholamines. The AR genes thus provide a paradigm for understanding the evolution of such structurally conserved yet functionally divergent families of receptor molecules.

Expansion of the alpha 2-adrenergic receptor family: cloning and characterization of a human alpha 2-adrenergic receptor subtype, the gene for which is located on chromosome 2.

Pharmacologic, biochemical, and genetic analyses have demonstrated the existence of multiple alpha 2-adrenergic receptor (alpha 2AR) subtypes. We have cloned a human alpha 2AR by using the polymerase chain reaction with oligonucleotide primers homologous to conserved regions of the previously cloned alpha 2ARs, the genes for which are located on human chromosomes 4 (C4) and 10 (C10). The deduced amino acid sequence encodes a protein of 450 amino acids whose putative topology is similar to that of the family of guanine nucleotide-binding protein-coupled receptors, but whose structure most closely resembles that of the alpha 2ARs. Competition curve analysis of the binding properties of the receptor expressed in COS-7 cells with a variety of adrenergic ligands demonstrates a unique alpha 2AR pharmacology. Hybridization with somatic cell hybrids shows that the gene for this receptor is located on chromosome 2. Northern blot analysis of various rat tissues shows expression in liver and kidney. The unique pharmacology and tissue localization of this receptor suggest that this is an alpha 2AR subtype not previously identified by classical pharmacological or ligand binding approaches.

Cloning of the cDNA for the human beta 1-adrenergic receptor.

Screening of a human placenta lambda gt11 library has led to the isolation of the cDNA for the human beta 1-adrenergic receptor (beta 1AR). Used as the probe was the human genomic clone termed G-21. This clone, which contains an intronless gene for a putative receptor, was previously isolated by virtue of its cross hybridization with the human beta 2-adrenergic receptor (beta 2AR). The 2.4-kilobase cDNA for the human beta 1AR encodes a protein of 477 amino acid residues that is 69% homologous with the avian beta AR but only 54% homologous with the human beta 2AR. This suggests that the avian gene encoding beta AR and the human gene encoding beta 1AR evolved from a common ancestral gene. RNA blot analysis indicates a message of 2.5 kilobases in rat tissues, with a pattern of tissue distribution consistent with beta 1AR binding. This pattern is quite distinct from the pattern obtained when the beta 2AR cDNA is used as a probe. Expression of receptor protein in Xenopus laevis oocytes conveys adenylate cyclase responsiveness to catecholamines with a typical beta 1AR specificity. This contrasts with the typical beta 2 subtype specificity observed when the human beta 2AR cDNA is expressed in this system. Mammalian beta 1AR and beta 2AR are thus products of distinct genes, both of which are apparently related to the putative G-21 receptor.

Sensitive and precise quantification of insulin-like mRNA expression in Caenorhabditis elegans.

Insulin-like signaling regulates developmental arrest, stress resistance and lifespan in the nematode Caenorhabditis elegans. However, the genome encodes 40 insulin-like peptides, and the regulation and function of individual peptides is largely uncharacterized. We used the nCounter platform to measure mRNA expression of all 40 insulin-like peptides as well as the insulin-like receptor daf-2, its transcriptional effector daf-16, and the daf-16 target gene sod-3. We validated the platform using 53 RNA samples previously characterized by high density oligonucleotide microarray analysis. For this set of genes and the standard nCounter protocol, sensitivity and precision were comparable between the two platforms. We optimized conditions of the nCounter assay by varying the mass of total RNA used for hybridization, thereby increasing sensitivity up to 50-fold and reducing the median coefficient of variation as much as 4-fold. We used deletion mutants to demonstrate specificity of the assay, and we used optimized conditions to assay insulin-like gene expression throughout the C. elegans life cycle. We detected expression for nearly all insulin-like genes and find that they are expressed in a variety of distinct patterns suggesting complexity of regulation and specificity of function. We identified insulin-like genes that are specifically expressed during developmental arrest, larval development, adulthood and embryogenesis. These results demonstrate that the nCounter platform provides a powerful approach to analyzing insulin-like gene expression dynamics, and they suggest hypotheses about the function of individual insulin-like genes.

Delivery and Scavenging of Nucleic Acids by Polycationic Polymers

Electrostatic interaction is a strong force that attracts positively and negatively charged molecules to each other. Such an interaction is formed between positively charged polycationic polymers and negatively charged nucleic acids. In this dissertation, the electrostatic attraction between polycationic polymers and nucleic acids is exploited for applications in oral gene delivery and nucleic acid scavenging. An enhanced nanoparticle for oral gene delivery of a human Factor IX (hFIX) plasmid is developed using the polycationic polysaccharide, chitosan (Ch), in combination with protamine sulfate (PS) to treat hemophilia B. For nucleic acid scavenging purposes, the development of an effective nucleic acid scavenging nanofiber platform is described for dampening hyper-inflammation and reducing the formation of biofilms.

Non-viral gene therapy may be an attractive alternative to chronic protein replacement therapy. Orally administered non-viral gene vectors have been investigated for more than one decade with little progress made beyond the initial studies. Oral administration has many benefits over intravenous injection including patient compliance and overall cost; however, effective oral gene delivery systems remain elusive. To date, only chitosan carriers have demonstrated successful oral gene delivery due to chitosan’s stability via the oral route. In this study, we increase the transfection efficiency of the chitosan gene carrier by adding protamine sulfate to the nanoparticle formulation. The addition of protamine sulfate to the chitosan nanoparticles results in up to 42x higher in vitro transfection efficiency than chitosan nanoparticles without protamine sulfate. Therapeutic levels of hFIX protein are detected after oral delivery of Ch/PS/phFIX nanoparticles in 5/12 mice in vivo, ranging from 3 -132 ng/mL, as compared to levels below 4 ng/mL in 1/12 mice given Ch/phFIX nanoparticles. These results indicate the protamine sulfate enhances the transfection efficiency of chitosan and should be considered as an effective ternary component for applications in oral gene delivery.

Dying cells release nucleic acids (NA) and NA-complexes that activate the inflammatory pathways of immune cells. Sustained activation of these pathways contributes to chronic inflammation related to autoimmune diseases including systemic lupus erythematosus, rheumatoid arthritis, and inflammatory bowel disease. Studies have shown that certain soluble, cationic polymers can scavenge extracellular nucleic acids and inhibit RNA-and DNA-mediated activation of Toll-like receptors (TLRs) and inflammation. In this study, the cationic polymers are incorporated onto insoluble nanofibers, enabling local scavenging of negatively charged pro-inflammatory species such as damage-associated molecular pattern (DAMP) molecules in the extracellular space, reducing cytotoxicity related to unwanted internalization of soluble cationic polymers. In vitro data show that electrospun nanofibers grafted with cationic polymers, termed nucleic acid scavenging nanofibers (NASFs), can scavenge nucleic acid-based agonists of TLR 3 and TLR 9 directly from serum and prevent the production of NF-ĸB, an immune system activating transcription factor while also demonstrating low cytotoxicity. NASFs formed from poly (styrene-alt-maleic anhydride) conjugated with 1.8 kDa branched polyethylenimine (bPEI) resulted in randomly aligned fibers with diameters of 486±9 nm. NASFs effectively eliminate the immune stimulating response of NA based agonists CpG (TLR 9) and poly (I:C) (TLR 3) while not affecting the activation caused by the non-nucleic acid TLR agonist pam3CSK4. Results in a more biologically relevant context of doxorubicin-induced cell death in RAW cells demonstrates that NASFs block ~25-40% of NF-ĸβ response in Ramos-Blue cells treated with RAW extracellular debris, ie DAMPs, following doxorubicin treatment. Together, these data demonstrate that the formation of cationic NASFs by a simple, replicable, modular technique is effective and that such NASFs are capable of modulating localized inflammatory responses.

An understandable way to clinically apply the NASF is as a wound bandage. Chronic wounds are a serious clinical problem that is attributed to an extended period of inflammation as well as the presence of biofilms. An NASF bandage can potentially have two benefits in the treatment of chronic wounds by reducing the inflammation and preventing biofilm formation. NASF can prevent biofilm formation by reducing the NA present in the wound bed, therefore removing large components of what the bacteria use to develop their biofilm matrix, the extracellular polymeric substance, without which the biofilm cannot develop. The NASF described above is used to show the effect of the nucleic acid scavenging technology on in vitro and in vivo biofilm formation of P. aeruginosa, S. aureus, and S. epidermidis biofilms. The in vitro studies demonstrated that the NASFs were able to significantly reduce the biofilm formation in all three bacterial strains. In vivo studies of the NASF on mouse wounds infected with biofilm show that the NASF retain their functionality and are able to scavenge DNA, RNA, and protein from the wound bed. The NASF remove DNA that are maintaining the inflammatory state of the open wound and contributing to the extracellular polymeric substance (EPS), such as mtDNA, and also removing proteins that are required for bacteria/biofilm formation and maintenance such as chaperonin, ribosomal proteins, succinyl CoA-ligase, and polymerases. However, the NASF are not successful at decreasing the wound healing time because their repeated application and removal disrupts the wound bed and removes proteins required for wound healing such as fibronectin, vibronectin, keratin, and plasminogen. Further optimization of NASF treatment duration and potential combination treatments should be tested to reduce the unwanted side effects of increased wound healing time.

Chiropteran types I and II interferon genes inferred from genome sequencing traces by a statistical gene-family assembler.

BACKGROUND: The rate of emergence of human pathogens is steadily increasing; most of these novel agents originate in wildlife. Bats, remarkably, are the natural reservoirs of many of the most pathogenic viruses in humans. There are two bat genome projects currently underway, a circumstance that promises to speed the discovery host factors important in the coevolution of bats with their viruses. These genomes, however, are not yet assembled and one of them will provide only low coverage, making the inference of most genes of immunological interest error-prone. Many more wildlife genome projects are underway and intend to provide only shallow coverage. RESULTS: We have developed a statistical method for the assembly of gene families from partial genomes. The method takes full advantage of the quality scores generated by base-calling software, incorporating them into a complete probabilistic error model, to overcome the limitation inherent in the inference of gene family members from partial sequence information. We validated the method by inferring the human IFNA genes from the genome trace archives, and used it to infer 61 type-I interferon genes, and single type-II interferon genes in the bats Pteropus vampyrus and Myotis lucifugus. We confirmed our inferences by direct cloning and sequencing of IFNA, IFNB, IFND, and IFNK in P. vampyrus, and by demonstrating transcription of some of the inferred genes by known interferon-inducing stimuli. CONCLUSION: The statistical trace assembler described here provides a reliable method for extracting information from the many available and forthcoming partial or shallow genome sequencing projects, thereby facilitating the study of a wider variety of organisms with ecological and biomedical significance to humans than would otherwise be possible.

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing.

BACKGROUND: There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important. RESULTS: Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage. CONCLUSIONS: We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.

Development of universal antidotes to control aptamer activity.

With an ever increasing number of people taking numerous medications, the need to safely administer drugs and limit unintended side effects has never been greater. Antidote control remains the most direct means to counteract acute side effects of drugs, but, unfortunately, it has been challenging and cost prohibitive to generate antidotes for most therapeutic agents. Here we describe the development of a set of antidote molecules that are capable of counteracting the effects of an entire class of therapeutic agents based upon aptamers. These universal antidotes exploit the fact that, when systemically administered, aptamers are the only free extracellular oligonucleotides found in circulation. We show that protein- and polymer-based molecules that capture oligonucleotides can reverse the activity of several aptamers in vitro and counteract aptamer activity in vivo. The availability of universal antidotes to control the activity of any aptamer suggests that aptamers may be a particularly safe class of therapeutics.

Mitochondrial mutations in adenoid cystic carcinoma of the salivary glands.

BACKGROUND: The MitoChip v2.0 resequencing array is an array-based technique allowing for accurate and complete sequencing of the mitochondrial genome. No studies have investigated mitochondrial mutation in salivary gland adenoid cystic carcinomas. METHODOLOGY: The entire mitochondrial genome of 22 salivary gland adenoid cystic carcinomas (ACC) of salivary glands and matched leukocyte DNA was sequenced to determine the frequency and distribution of mitochondrial mutations in ACC tumors. PRINCIPAL FINDINGS: Seventeen of 22 ACCs (77%) carried mitochondrial mutations, ranging in number from 1 to 37 mutations. A disproportionate number of mutations occurred in the D-loop. Twelve of 17 tumors (70.6%) carried mutations resulting in amino acid changes of translated proteins. Nine of 17 tumors (52.9%) with a mutation carried an amino acid changing mutation in the nicotinamide adenine dinucleotide dehydrogenase (NADH) complex. CONCLUSIONS/SIGNIFICANCE: Mitochondrial mutation is frequent in salivary ACCs. The high incidence of amino acid changing mutations implicates alterations in aerobic respiration in ACC carcinogenesis. D-loop mutations are of unclear significance, but may be associated with alterations in transcription or replication.

An Exploration into Fern Genome Space.

Ferns are one of the few remaining major clades of land plants for which a complete genome sequence is lacking. Knowledge of genome space in ferns will enable broad-scale comparative analyses of land plant genes and genomes, provide insights into genome evolution across green plants, and shed light on genetic and genomic features that characterize ferns, such as their high chromosome numbers and large genome sizes. As part of an initial exploration into fern genome space, we used a whole genome shotgun sequencing approach to obtain low-density coverage (∼0.4X to 2X) for six fern species from the Polypodiales (Ceratopteris, Pteridium, Polypodium, Cystopteris), Cyatheales (Plagiogyria), and Gleicheniales (Dipteris). We explore these data to characterize the proportion of the nuclear genome represented by repetitive sequences (including DNA transposons, retrotransposons, ribosomal DNA, and simple repeats) and protein-coding genes, and to extract chloroplast and mitochondrial genome sequences. Such initial sweeps of fern genomes can provide information useful for selecting a promising candidate fern species for whole genome sequencing. We also describe variation of genomic traits across our sample and highlight some differences and similarities in repeat structure between ferns and seed plants.

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey.

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson-Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1'-C1' distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1'-C1' distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5'-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.

Computational Methods for RNA Structure Validation and Improvement.

With increasing recognition of the roles RNA molecules and RNA/protein complexes play in an unexpected variety of biological processes, understanding of RNA structure-function relationships is of high current importance. To make clean biological interpretations from three-dimensional structures, it is imperative to have high-quality, accurate RNA crystal structures available, and the community has thoroughly embraced that goal. However, due to the many degrees of freedom inherent in RNA structure (especially for the backbone), it is a significant challenge to succeed in building accurate experimental models for RNA structures. This chapter describes the tools and techniques our research group and our collaborators have developed over the years to help RNA structural biologists both evaluate and achieve better accuracy. Expert analysis of large, high-resolution, quality-conscious RNA datasets provides the fundamental information that enables automated methods for robust and efficient error diagnosis in validating RNA structures at all resolutions. The even more crucial goal of correcting the diagnosed outliers has steadily developed toward highly effective, computationally based techniques. Automation enables solving complex issues in large RNA structures, but cannot circumvent the need for thoughtful examination of local details, and so we also provide some guidance for interpreting and acting on the results of current structure validation for RNA.